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SECTION  I 
INTRODUCTION 


One  of  the  more  popular  approaches  being  employed  in  rocket  nozzle 
structures  is  known  as  rosette  construction.  In  this  construction,  layers  of 
woven  fabric  composite  materials  are  wrapped  to  follow  spiral  trajectories 
in  the  cross-section  of  an  exit  cone  while  the  warp  or  fill  direction  generates 
a helix  about  the  longitudinal  axis  of  the  structure.  The  purpose  of  rosette 
construction  is  to  impart  a high  stiffness  in  the  radial  (thickness)  direction 
as  a means  of  improving  load- carrying  capacity  as  well  as  erosion  resistance 
as  compared  to  tangential  wrapping  of  the  layers. 

The  concept  of  rosette  construction  was  initially  devised  for  use  in  con- 
nection with  ablative  plastic  liners,  where  structural  performance  was  of 
minor  interest.  Consequently,  no  disciplined  efforts  were  undertaken  to 
relate  the  structural  response  to  the  established  effective  moduli  of  the  basic 
sheet  material.  The  concept  is  now  being  applied  in  the  construction  of  rocket 
nozzle  exit  cones,  where  the  structural  response  characteristics  are  critically 
important. 

In  this  work,  we  shall  develop  the  treatment  defining  the  exact  relations 
that  exist  among  the  various  geometric  parameters  in  a rosette  cylinder. 

We  shall  also  formulate  the  exact  differential  equations  (in  the  sense  of 
classical  anisotropic  elasticity  theory)  which  define  the  mechanical  response 
of  such  bodies  under  typical  loading  conditions  imposed  during  an  experimental 
characterization  program.  The  governing  equations  will  be  solved  numerically 
and  examples  presented  which  highlight  the  errors  involved  in  the  use  of  crude 
analytical  techniques.  It  will  be  shown  that,  in  general,  21  elastic  coefficients 
are  necessary  to  define  the  response  of  a rosette  cylinder.  Furthermore, 
the  21  coefficients  are  functions  of  position  in  the  medium.  While  these 
elastic  coefficients  are  not  mutually  independent  since  they  are  related  to  the 
stiffness  matrix  of  the  basic  sheet  material,  they  must  all  be  included  in  the 
theory  to  define  rigorous  solutions  of  boundary  value  problems.  Indeed, 
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even  the  proper  interpretation  of  "simple"  experiments,  such  as  axial  load- 
ing, torsion,  and  internal  pressure  depends  upon  such  solutions.  As  such 
experiments  are  fundamental  with  respect  to  the  structural  design  procedure, 
they  define  the  scope  of  the  present  work.  Included  in  the  treatment  will  be 
the  influence  of  uniform  "expansional"  strains  [1]  which  allow  us  to  consider 
the  stress  fields  induced  by  uniform  temperature  changes,  fabrication  curing 
conditions,  and  absorption  of  swelling  agents.  We  might  also  suggest  that  this 
work  is  basic  to  the  problem  of  optimization  of  the  material  and  geometric 
parameters  in  structural  elements,  such  as  rosette  cylinders  or  conical 
shells . 

SECTION  II 

GEOMETRIC  RELATIONS 

The  rosette  pattern  can  be  generated  by  starting  with  the  basic  sheet 
material  in  which  the  axes  of  elastic  symmetry  are  denoted  by  x^  (i  = 1,2,  3) 
as  shown  in  Figure  1(a).  The  sheet  is  then  distorted  to  follow  a certain 
cylindrical  spiral  trajectory,  defined  by  C and  C'  , while  the  deformed  axis 
Xj  assumes  the  helical  path  shown  in  Figure  1(b).  The  next  sheet,  initially 
oriented  as  in  Figure  1(a),  is  then  laid  adjacent  to  the  first  and  the  process 
is  continued  until  the  entire  area  A between  radii  r^  and  r^  is  completely 
filled.  As  the  specific  equation  of  curve  C is  not  known  apriori,  we  must 
define  the  pertinent  relations  existing  among  the  arc  length,  ply  thickness, 
and  inner  and  outer  radii  to  establish  this  unique  relation,  consistent  with 
the  continuity  requirement  (complete  filling  of  cross-sectional  area  A with 
no  gaps).  Finally,  in  order  to  develop  the  appropriate  constitutive  relations, 
the  path  of  the  original  x^  axis  within  the  rosette  cylinder  must  be  defined. 

We  begin  by  treating  the  trajectory  of  curve  C.  Consider  two  arbitrary 
spirals  C and  c'  which  intersect  a circle  of  radius  r at  points  A and  D as 
shown  in  Figure  2.  For  conceptual  purposes,  we  may  let  the  two  spirals 
represent  the  edges  of  a single  layer  of  the  basic  sheet  material,  which  we 
assume  is  completely  flexible.  Thus,  the  perpendicular  distance  between 
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the  two  spirals,  t,  must  be  a constant.  I.et  the  local  angles  between  the 
tangents  to  the  spirals  at  A and  B and  the  intersecting  circles  centered  at 
O be  denoted  by  q and  a'  , respectively.  We  shall  henceforth  refer  to  a as 
the  arc  angle.  Our  convention  is  that  the  angle  q is  measured  in  the  clock- 
wise direction  from  the  local  positive  9 direction  to  the  spiral  tangent.  By 
applying  the  law  of  sines  to  the  triangle  OAB,  we  get 

r sin  q = r ' sin  a1  ( 1) 

Since  there  is  not  a unique  trajectory  of  curve  C satisfying  the  continuity 
requirement,  we  select  the  class  of  curves  commonly  employed  in  rosette 
construction,  namely,  the  one  in  which  all  curves  intersect  a circle  of  given 
radius  at  the  same  arc  angle.  Therefore,  at  radial  distance  r'  , the  arc 
angle  on  C is  equal  to  q'  . Hence,  the  parametric  form  of  the  equation  * 
of  curve  C is  given  by  ( 1) , or 

r sin  a = r sin  a = const.  (2) 

o o 1 

where  a represents  the  arc  angle  corresponding  to  the  point  (r  ,9  ) in 
o o o 

polar  coordinates. 

To  continue  our  study  of  the  spiral  configuration,  we  consider  the 
differential  geometric  relation  (see  Figure  3) 

dr  = - r tan  a d9  (3) 

Eliminating  the  dependence  on  r in  favor  of  a through  (2)  and  (3),  and  inte- 
grating yields 

9-0  = a + cot  a - q - cot  a (4) 

o o o ' ’ 

And  finally,  using  (2)  again,  we  can  express  the  equation  of  C in  terms  of 

polar  coordinates  as 

+ 

1 This  equation,  as  well  as  several  others  presented  in  this  section  have 
been  derived  by  Mamrol  [2]  using  a different  approach. 
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r -r  sm  a ) 

V o of 


r sm  a 
o o 


where  the  ambiguous  sign  is  negative  if  0<ao<  ff/2  and  positive  if  n/ 2<a^<  it. 

In  order  to  fabricate  the  rosette,  we  need  to  compute  the  arc  length  of 
the  spiral.  From  Figure  3,  we  observe  that  infinitesimal  arc  length  ds  is 
given  by 


Substituting  for  dr  from  (2)  into  (6)  and  integrating,  we  get 


2sina 


. 2 
sm  a 

c 

. 2 
sm  a 


r 2 

S = TT-2 ( 1 - £o  ) 

2sina  2 

o r 

o 

both  of  which  represent  the  arc  length  from  (r0>®0)  to  an  arbitrary  point 


(r,9), 


In  order  to  define  the  relation  between  two  arbitrary  spirals,  we  again 


return  to  Figure  2,  noting  that 

t sin  a'  = r sin  (a-a'  ) 

by  use  of  the  law  of  sines.  Solving  (9)  for  a'  gives 


cot  a'  = cot  a + 


rsina 


By  applying  (4)  to  points  D and  B and  taking  the  difference,  we  find  that 

^ = a'  + cot  a'  - a - cot  a 

Thus  the  central  angle  <J>  subtended  by  arc  AD  is  given  by 


rsina  r sina 
o o 


on  two  given  spirals,  the  points  being  at  the  same  radial  distance  from  O,  is 
, constant.  Therefore,  given  the  trajectory  of  one  spiral  within  the  rosette 

[ pattern,  any  other  spiral  is  generated  by  a rigid  body  rotation  about  the  longi- 

| tudinal  axis  through  point  O.  The  magnitude  of  the  angle  of  rotation  is  given 

\ by  (12),  where  t is  the  perpendicular  distance  separating  the  spirals.  Hence, 

the  number  of  layers  N within  the  cylinder  is 

2 n r s in  q 

N = (13) 

1 

where  t is  the  thickness  of  a single  layer.  Satisfaction  of  Eqs.  (13)  and  (4) 
or  (5)  implies  satisfaction  of  the  continuity  requirement.  Finally,  given  the 
equation  of  one  spiral,  i.e.,  (4),  the  equation  of  any  arbitrary  spiral  can  be 
expressed  as 

9-0  - <t>  = a + cot  a - q - cot  a (14) 

o o o 

where  <J> depends  on  dimension  t through  Eq.  (12). 


In  order  to  relate  the  helical  angle  o>  to  the  respective  angle  in  the 
basic  sheet,  <t>  , we  refer  to  Figure  4,  where  the  basic  sheet  is  illustrated  in 
4(a)  and  the  top  view  of  the  sheet  in  the  rosette  cylinder  is  shown  in  4(b). 
First,  we  observe  the  relation 


z = £ cot  <p 

where  £ can  be  evaluated  from  (8).  Thus,  we  get 


2r  sin  a 
o o 


Now,  helical  angle  co  is  defined  by 


(15) 


(16) 


COS  (O 


dz 

dR 


(17) 
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where  dR  is  the  infinitesimal  arc  length  along  the  deformed  line  QP  in  the 
rosette  cylinder,  so  that 


dR 


, 2 . 2 

= dz  + dr 


2 , 2 
r d0 


(18) 


Taking  differentials  of  (4)  and  ( 1 6) , using  (2),  and  substituting  the  results  into 
( 1 8)  gives 


dR 


- rdr 


r sin  a sincf) 
o o 


(19) 


Use  of  (17)  and  (19),  in  conjunction  with  the  evaluation  of  dz  from  (16)  yields 

cos  (D  = cos  4>  (20) 

Hence,  in  the  transformation  of  the  basic  sheet  into  the  rosette  pattern,  angles 
with  respect  to  the  longitudinal  direction  are  preserved.  This  implies  that 
the  angle  between  line  segments  in  the  basic  sheet  intersecting  at  an  arbitrary 
point  F is  carried  into  the  identical  angle  between  the  corresponding  tangents 
at  F'  through  the  present  transformation,  where  F'  is  the  transformed  posi- 
tion of  F . 


SECTION  III 
STRESS  ANALYSIS 

In  this  section,  we  shall  study  the  elastic  response  of  a rosette  cylinder 
under  the  application  of  an  axisymmetric  loading  system,  such  as  that  usually 
imposed  in  experimental  characterization,  in  addition  to  uniform  expansional 
strains  [lj.  The  latter  permit  the  treatment  of  thermal  stresses  induced  by 
uniform  temperature  changes,  curing  stresses  [3,4],  and  swelling  stresses 
caused  by  absorption  of  fluid.  We  assume  that  the  stress  field  is  independent 
of  the  axial  coordinate  z.  Thus,  the  localized  disturbances  near  the  ends  of 
the  specimen  under  imposed  displacement  boundary  conditions  are  excluded  in 
the  present  work.  As  mentioned  earlier,  the  rosette  cylinder  is  constructed 
of  numerous  layers  of  a composite  material.  In  what  follows,  each  layer  is 
treated  as  a homogeneous,  anisotropic  material  characterized  by  its  effective 
moduli  [5].  Thus  the  mechanical  properties  which  govern  the  response  are 
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the  effective  stiffness  coefficients  and  expansional  strains  of  the  basic  sheet 
material.  In  practical  application..,  the  basic  sheet  material  is  orthotropic 
(9  elastic  moduli  and  3 expansional  strains),  so  we  shall  consider  this  mater- 
ial class  in  the  present  work,  although  treatment  of  materials  possessing 
arbitrary  anisotropy  can  be  incorporated  without  difficulty. 


It  is  important  to  recall  that  the  rosette  structure  involves  a system  of 
layers,  the  trajectories  of  which  are  mutually  related  via  rigid  body  rotations. 
Furthermore,  the  helical  angles  defined  by  the  various  spiral  layers  are  all 
identical.  These  facts  imply  that  the  material  structure  is  axisymmetric,  i.e., 
while  the  moduli  are  functions  of  radial  position,  they  are  independent  of  &. 

As  the  loading  is  also  axisymmetric,  the  most  convenient  coordinate  system 
to  employ  in  connection  with  the  stress  analysis  is  that  of  cylindrical  coordin- 
ates r,  0,  z,  rather  than  coordinates  defined  by  the  spiral  trajectories. 


The  stress  components  in  the  present  class  of  boundary  value  problems  i 

1 

are  functions  of  radial  position  only.  It  follows  that  the  strain  distribution  j 

will  only  be  a function  of  r.  By  use  of  the  strain- displacement  relations  of 
linear  elasticity,  i.e.. 


ee  = r(U+'V  • % = W’z 


(21) 


y „ = v,  + — (u,  -v)  , y = u,  + w, 

'r9  r r 0 ' rrz  z r 


y n = v,  + - w 
z9  z r 9 


where  u,  v,  and  w represent  the  r,9,  and  z components  of  displacement, 

respectively,  and  commas  denote  differentiation,  it  can  be  shown  that  the 

2 

general  form  of  the  displacement  field  is  given  by 

u = U(r) 

v = Arz  + V(r)  (22) 

w = ez  + W(r) 


We  assume  that  the  displacement  field  is  continuous, 
slits  in  the  cylinder. 


i.e., 


there  are  no 
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where  € and  A are  constants.  By  substituting  (22)  into  (21),  the  strain 
field  can  be  expressed  as 


e = U,  , en  = U/r  , e = e 
r r 9 z 


Y a = V,  -V/r  , y = W,  , Y = Ai- 
rS  r rz  r zG 


(23) 


The  stress  equations  of  equilibrium  for  axisymmetric  problems  take  the 
form 

<r  + — (cr  -<r  ) = 0 

r,  r r r 9 


, , 

rz,  r r rz 


(24) 


r + — r =0 
r9 , 9 r r9 


Rewriting  the  first  of  (24)  and  integrating  the  last  two  directly  yields 


(rcr  ) = <rn 

r , r 9 


r r = const.  = B (25) 

rz 

2 

r r „ = const.  ^ D 

rG 

The  governing  field  equations  are  thus  defined  by  (23),  (24),  and  the  constitutive 
equations,  i.  e. , 


cr. . 


^ijkf^ki  6kf^ 


(26) 


where  are  the  components  of  the  mathematical  strain  tensor,  in  contrast 
to  those  of  (231  in  which  the  shear  strain  components  are  engineering  strains, 
and  e^  are  the  components  of  the  mathematical  expansional  strain  tensor. 

In  order  to  develop  the  proper  form  of  the  rosette  stiffness  coefficients, 
we  refer  to  the  orthogonal  coordinate  systems  shown  in  Figure  5.  Coordinates 
x.(i  = l,2,3)  are  oriented  along  the  axes  of  elastic  symmetry  of  the  basic  sheet 


8 


1 


material,  x.  are  tangent  and  normal  to  an  arbitrary  spiral  path,  and  x.  are 

the  cylindrical  coordinate  axes  9,r,z,  where  each  coordinate  system  is  right- 

handed.  Letting  x^  and  y^  represent  two  arbitrary  orthogonal  coordinate 

systems  and  C..,  and  B. , the  respective  stiffness  coefficients  in  the  two 
ijkf  ljkf 

systems,  the  transformation  equations  are  given  by 


B a .a  .a  a C 

ijkf  = pi  qj  rk  sf  pqrs 


(27) 


where  a.,  is  the  cosine  of  the  angle  between  x.  and  y..  The  stiffness  coeffi- 
*J  i J 

cients  in  the  x!  system,  C!  , may  thus  be  computed  by  applying  (27)  for  the 
i ljkf 

transformation  from  x.  to  x.,  followed  by  that  from  x,  to  x1.  The  results  of 

ii  ii 

these  transformations  are  given  in  the  appendix.  It  will  suffice  at  this  point 
to  state  that,  in  general,  all  elements  of  C\  are  non- zero,  but  they  may 

all  be  computed  in  terms  of  the  (nine)  stiffness  coefficients  of  the  basic  sheet 
Finally,  to  save  writing,  we  shall  subsequently  employ  the 
standard  contracted  notation  [6]  to  represent  the  stiffness  coefficients  in 


material,  C.., 
ljkf 


terms  of  double  indices,  e.g.,  C . Thus,  the  9 independent  components  C.. 
transform  into  13  non- zero  components  C..,  which  in  turn  are  carried  into 
21  components  C! ..  The  constitutive  relations  now  assume  the  form 


c r 

C C'  C'  C'  C'C'  C'  , 

r A 

X 

e 

11  12  13  14  15  16 

U/r  - e9 

\ 

\ 

0* 

C'  C!  C'  C'  C'  , 

: 

r 

22  23  24  25  26 

U - e 

5 

, r r 

-) 

<r 

C'  C'  C'  C'  . 

e - e 

J 

z 

33  34  35  36 

z 

(28)  j 

r 

SYMM.  C'  C'  C' 

W - e 

! 

rz 

44  45  46 

, r rz 

rz0 

C'  C'  , 
55  56 

Ar  - e „ 
z9 

! 

r A 

C' 

V - V/r-e 

l rel 

[ 66 

, r r9 

\ 
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where 


C!.<r) 


The  expansional  strain  components,  denoted  by  the  symbol  e,  are  also  func- 
tions of  r (except  for  e , which  is  a constant)  and  are  related  to  the  sheet 

z 

expansional  strains  e.(i=  1,2,3)  through  the  standard  transformation  of 
(engineering)  strain  equations  (see  appendix).  For  the  case  of  thermal 
expansion,  the  latter  strains  are  given  by 

e.  = a.  T (30) 

i i 

where  a.  are  the  coefficients  of  linear  thermal  expansion  and  T is  the  differ- 
ence between  ambient  temperature  and  that  in  the  stress- free  state.  For 
fluid  absorption,  e.  are  defined  by 

e.  = p.m  (31) 

l i 

where  p.  are  the  free  swelling  strain  components  per  unit  volume  of  fluid  in 
the  medium,  while  m is  the  total  volume  of  fluid  absorbed. 

Upon  evaluating  the  stiffness  coefficients  C\.  and  substituting  (28)  into 
(25),  we  obtain  3 ordinary  differential  equations  in  terms  of  dependent  varia- 
bles U,  V,  and  W.  Introducing  the  expression 

r sin  a 
o o 


these  governing  field  equations  become 

G r U + G U + G.  U/r  + G r V +G,(V  -V/r)+Gr\V 
3 ,rr  2 ,r  1 7 , rr  6 , r 5 , rr 


+ G . W 
4 , r 


G„  + G e + G Ar 
8 o 9 


Ci4U.r*C'l4U/r+Ci6,V.r-V,r,+  Ci4  W,  r ' H , 4B/ r 4C  3 5 p.  * APr(l-P 

i 2 1/22 

Ck6U.r+C16U/r4Ci6<V.r-V/r)tCi6W.r  <C8+  °o*,|SI  > -c'56AriWr 


in  the  region  r <r<r  . The  various  coefficients  in  (33)  are,  in  general,  func- 
1 o 

tions  of  r and  are  given  in  the  appendix. 
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The  boundary  conditions  required  to  complete  the  definition  of  the 
present  boundary  value  problem  are  given  by 


cr  (r  ) = cr 

r o o 


V(a)  = V 


%<ri>  = ri 


W(b)  = W 


where  cr  , cr  , V , and  W are  prescribed  constants,  and  a and  b are  arbi- 
o 1 o o 

trary  values  of  r within  the  medium.  The  constants  V and  W are  simply 

o o 

required  to  define  the  rigid  body  displacement  components.  Finally,  the 

constants  e.  A,  B,  and  D must  be  prescribed.  Note  that  B and  D completely 

define  shear  stresses  r and  r according  to  (25)  and  that  the  latter 

rz  r9 

normally  vanish  under  typical  experimental  loading  conditions. 

The  remaining  quantities  of  interest  are  the  axial  force  P and  torque 

o 


which  can  be  expressed  as 


P = 2 
o 


’ / 


(j  rdr 
z 


" / 


2, 

r r dr 
z6 


The  method  of  solution  adopted  here  requires  the  prescription  of  e and  A, 

rather  than  applied  loads  P and  T . However,  solutions  for  specified 

- o o 

values  of  P^  and  T^  can  be  developed  through  a simple  influence  function 
approach. 

SECTION  IV 

SOLUTION  OF  THE  BOUNDARY  VALUE  PROBLEM 


The  complexity  of  the  coefficients  involved  in  field  equations  (33),  in 
particular,  the  appearance  of  the  irrational  functions,  precludes  efficient 
application  of  the  standard  series  methods  of  solution  for  ordinary  differential 
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equations.  Consequently,  we  have  adopted  the  finite-difference  approach, 
incorporating  the  use  of  higher  order  difference  representations.  Forward 
differences  have  been  utilized  throughout  the  solution  with  a maximum  of 
4 nodal  point  functions  being  employed  in  the  difference  approximation.  Since 
the  first  of  (33)  involves  second  derivatives  of  the  displacement  functions 
and  we  wish  to  satisfy  the  field  equations  at  all  points  in  the  region,  including 
the  boundaries,  we  have  expanded  the  region  by  adding  two  nodal  points  out- 
side the  right  hand  boundary  r = r . 

Consider  the  series  of  equi-spaced  nodal  points  x = 0,h,  2h,  3h,  or 

x.  = ih  (i=0,  1,2,3)  where  x = r+d  and  d is  a constant.  Letting  y represent 
i o 

an  arbitrary  function  of  x (or  r)  and  letting  y^  = y(x^),  we  approximate  the 
first  and  second  derivatives  of  y,  in  the  event  that  all  4 nodal  points  lie  in  the 
(expanded)  region,  by 

hyo  T yo  + 3yl‘  2y2  +Ty3  • 

r < r - h (36) 

2 ° 
h y„'  * 2y„  • 5yl  + 4y2  - y3  • 

where  primes  denote  derivatives.  Equations  (36)  are  easily  derived  by 

expanding  y , y , y in  Taylor  series  about  y truncating  each  series  after 
I 5 o, 

the  cubic  term  in  h,  and  then  solving  for  y'  and  y".  At  r = r , three  point 

o o o 

difference  approximations  are  employed 


h y;  35  - I yo  + 2yl  - 2 y2 

h2yn  * yo  - 2yt  + y2  , 


r = r 

o 


(37) 


and  finally,  at  r = r^  + h,  we  put 

h y'0  - yi  ■ yo  * r = r0  + h 


(38) 
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Thus,  letting  N represent  the  total  number  of  nodal  points  employed  to  sub- 
divide the  expanded  region,  we  write  the  finite-difference  approximation  of 
the  first  of  (33)  at  N-2  nodal  points  and  that  of  the  second  and  third  of  (33) 
at  N-l  points.  Applying  the  finite- difference  representations  of  the  traction 
boundary  conditions  in  (34)  after  substituting  the  expression  for  a given  by 
(28)  and  the  two  conditions  for  V and  W from  (34),  we  arrive  at  3N  equations 
to  be  solved  for  the  3N  unknowns,  i.e.,  U,  V,  and  W at  the  N nodal  points. 
Utilizing  a computer  routine  for  sparsely  populated  matrices,  the  solution  of 
a given  problem  is  executed  in  a matter  of  seconds  on  a CDC6600. 


Upon  solving  for  the  nodal  displacement  functions,  the  stress  and 
strain  fields  are  defined  through  the  use  of  eqs.  (28)  and  (23),  respectively, 
where  derivatives  are  evaluated  via  the  appropriate  relation  among  (36)-(38). 
Since  the  solution  technique  is  an  approximate  one,  the  degree  of  approxima- 
tion may  be  assessed  by  writing  a system  of  "check"  equations,  consisting  of 
(25),  the  compatibility  relation. 


e = e„  + r e„ 
r 6 6,  r 


(39) 


and  the  identity 


r c 
o o 


Vi 


which  represents  the  integral  of  the  first  of  (25). 

SECTION  V 

SAMPLE  PROBLEM 


(40) 


In  this  section  we  shall  consider  the  stress  field  with'Ln  a rosette 
cylinder  in  which  the  basic  sheet  material  is  a high  modulus  unidirectional 
graphite/epoxy  composite.  Letting  x^  represent  the  direction  parallel  to 
the  fibers,  we  assume  the  following  basic  sheet  moduli  (psi), 

Cn  = 20.2  x 106  , C = CJ3  = 4.72  x 105 


23 


= 3.84  x 10' 


22 


= C 


33 


= 1.50  x 10 


1 3 


(41) 


(41  cont.) 


C„  „ = 5.  60  x 10^  , C = C,.  = 7.0  x 105 

44  55  66 

The  cylinder  is  under  an  internal  pressure  = -1000  psi,  while  cr  = B = 

D = 0,  i.e.,  the  outer  surface  is  traction- free  and  the  shear  traction  on  the 

inner  surface  vanishes.  We  also  let  e - A = 0,  which  imply  that  axial  force 

P and  torque  T , in  general,  are  different  from  zero.  The  inner  and  outer 
o o 

radii  are  3"  and  4",  respectively,  and  the  expansional  strains  e.  are  assumed 
to  vanish. 

Three  specific  material  configurations  will  be  treated.  In  the  first 

(rosette),  we  put  a = JT / 1 2 and  4>  = ff/3.  Two  other  configurations  which 

may  be  employed  to  crudely  approximate  the  rosette  structure  have  also  been 

considered.  These  are  represented  by  the  parameters  a = 0 , <t>  - */3 

o 

(helical- wound  cylinder)  and  aQ  = 0,  <f>  = n/2  (hoop-wound  cylinder).  In  the 
latter  two  configurations,  exact  solutions  for  the  present  class  of  boundary 
value  problems  are  available  (see  ra  ).  Such  solutions  were  also  treated  by 
use  of  the  present  approach  to  serve  in  checking  the  computer  code.  In  the 
results  presented  here,  403  node  points  were  utilized  in  each  solution,  which 
is  sufficient  to  achieve  a precision  of  approximately  4 significant  figures  for 
the  computed  stress  components. 


Distributions  of  the  four  non- vanishing  stress  components,  cr^,  cr^,  r^g, 
and  <y  , throughout  the  domain  are  given  in  Figures  6-9  for  the  foregoing 
3 configurations.  Clearly,  the  helical- wound  and  hoop- wound  cylinder  models 
represent  very  poor  approximations  to  the  actual  response,  although  smaller 
errors  than  those  indicated  here  may  be  expected  in  systems  of  less  extreme 
anisotropy  than  the  present  case.  From  these  curves,  we  also  observe  that 
the  rosette  structure  effects  a considerable  smoothing  of  the  stress  field,  at 
least  under  the  conditions  treated  here,  i.e.,  stress  concentrations  in  cr  , 
cr  , and  r are  close  to  unity.  While  the  cylindrical  (r,9,z)  components  of 
the  stress  tensor  are  not  particularly  informative  with  regard  to  failure 
prediction  in  orthotropic  bodies,  the  stress  components  in  the  elastic  symme- 
try coordinates  x,  (see  Figure  5)  can  be  easily  determined  through  the 
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standard  stress  transformation  equations.  Consequently,  coupled  with  an 
appropriate  failure  criterion  for  the  basic  sheet  material,  the  present  solu- 
tion provides  a means  of  optimizing  the  material  and  geometric  parameters 
in  rosette  cylindrical  bodies  under  axisymmetric  loading. 

SECTION  VI 

CONCLUSIONS 

We  have  derived  the  pertinent  relations  which  describe  the  geometry  of 
the  spiral  paths  traversed  by  the  layers  within  a rosette  cylindrical  body.  It 
has  been  demonstrated  that  prescription  of  inner  and  outer  radii,  r^  and  r^, 
respectively,  arc  angle  a at  one  point,  and  layer  thickness  t^  completely 
defines  the  appropriate  geometric  pattern.  The  spiral  trajectories  of  the 
various  layer  interfaces  are  interrelated  through  simple  rigid  body  rotations, 
while  helical  angle  is  constant.  These  facts  imply  axisymmetry  of  the 
distribution  of  heterogeneous  effective  stiffness  coefficients  C!,. 

Taking  advantage  of  the  axisymmetric  material  structure,  we  have 
formulated  the  governing  equations  of  elasticity  for  the  rosette  cylinder  under 
boundary  conditions  and  expans ional  strains  which  are  independent  of  9. 
Although  the  stiffness  coefficient  matrix  C.’.  is  fully  populated  (21  constants) 
in  general,  these  coefficients  only  depend  upon  the  (orthotropic)  moduli  of 
the  basic  sheet  material  and  the  angles  a and  <f>  . The  finite- difference  method 
has  led  to  the  development  of  an  efficient  computer  program  for  the  general 
solution  of  the  axisymmetric  class  of  boundary  value  problems.  An  example 
solution  has  demonstrated  the  tremendous  smoothing  influence  of  the  rosette 
construction  on  the  elastic  stress  field,  which  will  be  an  important  feature  of 
this  concept  in  structural  applications. 
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APPENDIX 


The  components  of  the  stiffness  matrix  (see  Figure  5)  are  defined 


by  substituting 


a = a = a = a =0 
12  22  31  33 


into  (27).  Upon  switching  to  standard  contracted  notation,  we  get 

Si  =Chn4+C22m4+2<Ci2  + 2C66>MV 


C12'C,3n2+C23mZ 


C13  =(CU+C22-  4C66)M2N2+C12(M4+N4) 


C 1 5 = [CnN2-  C22M2t(C12+2C66)(M2-  N2)]  MN 


C22  = °33 


C23=C13m2+C23n2 


C25=(C13-C23>MN 


C33  = cuM4+C22N43  2(Ci2+2C66,M2N2 


C35  = [C11m2-  C22n2+(C12+2C66>(n2-  “'’J  MN 


c,,  = c m2+  c . .n2 

44  55  44 


C46  = (C55-  C44>MN 


c«  (C11,C22>mV-  2C12M2N2tC66,MZ-  nV 


C.  c N24C  „M2 
bb  S5  44 


Ci4  = Ci6  = ° 


(i  = 1,2,3,  5) 


Similarly,  the  transformation  from  x.  to  x^,  upon  substituting 

2 1/2 

a j j = a22  = cos  Q = ( 1 - p ) , a2i  = -al2  = sin  q e (3  , = 1 


a 1 3 r a23  = a31  = a32  ° 


into  (27)  yields 


c;,  = C,  ,( l-p2)2  +C,,p4  + 2(C12+ 2Cfe6)p2(l-p2) 


11  11 


C’l2  =IC1I,C22-  4Cfc6,P2,1-<52,+  C|2<2P4-  2p2  +1) 


C13  ° Ci3<‘-P2>+  C23P2 


Ci4  ' f(ZC46-  C15)(1-PZ)  -Cz/]p 

? — 2 21/2 

C’l5  =[Cl5"-fi  >+<C25+  2C46)P  I(1-P  > 

Ci6=l=n(pZ-1'+522p2+<h2+2566,(!-2|52,1|5<1-p2)‘/2 

Ci2  = S2(l-|32)2t  2(C12+  2C66)P2(1-P2) 


cb  = cnp2+c23<1-fj2) 


C24=  [-c15(52t(C25+2C46)(p2-ll]U 

_ — 2—  2 21/2 
C25  = «C15-  2C46)P  +C25  d-M)0-M 


— 2—  2 — — 2 21/2 
CL,  f-  cnP  + C22(i-p  )4<C1?+  2CAA)(2p  -1)]  (3(1- p ) 
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C 3 3 = C 3 3 


C34  - -C35P 


2 1/2 

C35^C3S(1-P  > 


c,36  = ,c23-cn’p"',i2|‘/2 


CU  =C44,1-|i2'  +C5/ 


C45  = (C44-  c55'PO-P2'1/2 


(A- 4 cont, ) 


— — 2 21/2 
Ci6  =[(C15-  C25l(!  *C46<‘-2P  11  (1-p  > 


c^5  = C44P2+  c55<’-PZ| 


C56  =f  (CZ5-  C15>"-P2>+  C46(I-2p2,  j p 


C66  = (C11  tC22-  2C12)p2U-p2)+C66(1-2(52,2 


where  we  have  assumed  that  q<  */2.  Note  that  C.,  in  the  above 
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expressions  are  all  constants. 

Substituting  (A-l)  and  (A-3)  into  the  standard  strain  transformation 
equations  leads  to 

2 2 2 2 2 
e = (e  - e N - e M )p  + e N + e_M 
9 3 12  12 

222 

e =(e,N  + e_M  - e )6  +e 
r 1 2 3^3 


2 2 

e = e,  M + e_N 

7,  1 2 


(A- 5) 
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( A-  5 cont. ) 


where  we  recall  that  e are  the  engineering  expansional  strain  components. 

y 

FinalLy,  the  remaining  functions  appearing  in  field  equations  (33)  are 
given  by 

co  = 5n-  S3 

Gi  = 3F/-  ZF2P2-en 


C22‘  Cll'  G1 


= V4.2F/t  c22 


(A-6) 


= (F,pScJ5-  ZC46)p 


= -(F3P2  + C25t2C46»! 


= ( -3F1p4t(F14  3F2ip2  +Cn-  C22J|3U-p2)‘1/2 


G. 


- (-FjP  + F2)p(i-p  ) 


2 1/2 


2 2 - - 


.2  .2 


(C12'  Cll)(eiN  +62M  ^C22-Cl2)e3+(C23-Cl3)(eiM  +e2N  > 


+2(C25-C15)(ere2)MN 


2 — — 2-1/2 
C9  = (-2F3P  +C15-2C25)(l-p  ) 
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(A-  6 cont. ) 


H *'  -fc15(e1N2+e2M2)+C25e3+C35(eiM2+e2N2)+2C55(ere2)MNj|3 


wh  ere 


f,  = cu  +c22-  2ci2-4c66 


F2  = C22-C12-2C66 


(A- 7) 


r3  = C15-  C25-  2C46 


This  completes  the  definition  of  all  functions  appearing  in  the  governing 
field  equations  (33). 
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